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Abstract. We prove convergence of a fully discrete finite difference scheme 
for the Korteweg— de Vries equation. Both the decaying case on the full line 
and the periodic case are considered. If the initial data u|t=o = uo is of high 
regularity, uq £ H 3 (R), the scheme is shown to converge to a classical solution, 
and if the regularity of the initial data is smaller, uo £ L 2 (K), then the scheme 
converges strongly in L 2 (0, T; L 2 oc (R)) to a weak solution. 



1. Introduction 

The Korteweg-de Vries (KdV among friends and foes) equation, which reads 
(1.1) u t + uu x + u xxx = 0, 

has been studied extensively since its first analysis in 1895 by Korteweg and de Vries. 
Apart from applications as a model for shallow water waves, the KdV equation 
has maintained a pivotal role in several branches of mathematics. We here focus 
on the derivation of convergent numerical methods for the initial value problem 



where the equation (1.1) is augmented by initial data u\ t= o = u . The problem of 
analyzing convergent numerical schemes is of course intimately connected with the 
mathematical properties of the Cauchy problem for the KdV equation, which has 
undergone a tremendous development the last two decades, see, e.g., |TT1 [14] and 
the references therein. We will not be able to discuss this literature here, but only 
refer to the parts that are pertinent to the current paper. 

In this paper we analyze the implicit finite difference scheme 

(1.2) v% +1 =tq-AtuqDu^-AtD 2 + D_u^ +x , n G No, j S Z, 

where u" ~ u(jAx, nAt), and Ax, At are small discretization parameters. Further- 
more, D and D± denote symmetric and forward/backward (spatial) finite differ- 
ences, respectively, and u denotes a spatial average. Two results are proven, both 
for the full line and the periodic case: (1) In the case of initial data uq € i/ 3 (R), 



we show (see Theorem 3.3 and Remark 3.4) that the approximation (1.2| converges 
uniformly as Ai -> with At = O (Aarjin C(R x [0, T]) for any positive T to the 
unique solution of the KdV equation. (2) When the initial data uq E L 2 (K), we 



prove that (see Theorem 4.3 ) that the approximation converges strongly as Ax — ¥ 
with At = (D (Ax 2 ) in LF(Q,T\ Lf oc (K)) to a weak solution of the KdV equation. 

An interesting fact, and rarely referred to in the current literature, is that the first 
mathematical proof of existence and uniqueness of solutions of the KdV equation, 
was accomplished by Sjoberg |16j in 1970, using a finite difference approximation 
very much in the spirit considered here. His proof is valid for initial data that are 
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periodic and with square integrable third derivative, that is, uq{x + 1) = uo(x) for 
x £ [0, 1] and u'q £ L 2 ([0, 1]). Sjoberg's uniqueness proof still is the standard one, 
using the Gronwall inequality. His approach is based on a semi-discrete approx- 
imation where one discretizes the spatial variable, thereby reducing the equation 
to a system of ordinary differential equations. However, we stress that for numer- 
ical computations also this set of ordinary differential equations will have to be 
discretized in order to be solved. Thus in order to have a completely satisfactory 
numerical method, one seeks a fully discrete scheme that reduces the actual com- 
putation to a solution of a finite set of algebraic equations. This is accomplished 
in the present paper, both in the periodic case and on the full line. 

There has been a number of papers involving the numerical computation of 
solutions of the Cauchy problem, starting with the landmark paper by Zabusky and 
Kruskal [20] , where they discovered the permanence of solitons (the term "solitons" 
being coined in the same paper) for the KdV equation using numerical techniques. 
However, we will here focus on papers that discuss numerical methods per se. 

A popular numerical approach has been the application of various spectral meth- 
ods. Little is known rigorously about the convergence of these methods. For a sur- 
vey and a comparison, see [15]. See also [12]. Multisymplectic schemes have been 
studied in [3] (see also references therein) . There exist convergence proofs for finite 
element methods for the KdV equation, see [HI 2J However, the resulting 
schemes tend to be quite different from finite difference schemes derived ab initio. 

The numerical computation of solutions of the KdV equation is rather capricious. 
Two competing effects are involved, namely the nonlinear convective term uu x , 
which in the context of the Burgers equation Ut + uu x = yields infinite gradients 
in finite time even for smooth data, and the linear dispersive term u xxxi which in 
the Airy equation u t + u xxx = produces hard-to-compute dispersive waves, and 
these two effects combined makes it difficult to obtain accurate and fast numerical 
methods. Indeed, any initial data for the Burgers equation that is decreasing in 
a small neighborhood, will develop infinite gradients in finite time, while the Airy 
equation preserves all Sobolov norms while creating many oscillatory waves. Most 
finite difference schemes will consist of a sum of two terms, one discretizing the 
convective term and one discretizing the dispersive term. These two effects will 
have to balance each other, as it is known that the KdV equation itself keeps the 
Sobolov norm H 8 bounded; from [6] we know that if uo £ H S (M), with s > 3, then 
the solution satisfies ||w(£)IIh=(r) < Ct,u f° r t G [O,? 1 ]. This dichotomy between 
these two effects is brought to the forefront in the method of operator splitting. 
Here the two equations, the Burgers equations and the Airy equation, are solved 
sequentially for a small time step. This procedure is iterated, and as the time 
step converges to zero, the approximation converges to the actual solution. In the 
KdV context operator splitting was introduced by Tappert [18], a Lax-Wendroff 
theorem was proved in [7J, and convergence of the operator splitting technique 
proved in [8J [TTJ [SJ [TU] . Our approach here is a finite difference method which can 
also be viewed as an operator splitting method. 

Recently, a semi-discrete scheme for the generalized KdV equation was shown 
to converge in L 2 oc for initial data in L 2 \\\. However, the scheme analyzed here, 
which in contrast to the scheme in [I], does not involve an explicit fourth order 
stabilizing term, and we show convergence for non-smooth initial data. 

The rest of this paper is organized as follows: In Section [2]we present the nec- 
essary notation and define the numerical scheme. In Section[3] we show the con- 
vergence of the scheme for initial data in H 3 (R) , while in Section [4] we show the 
convergence to a weak solution if the initial data is in L 2 (M). In Section [5] we 
exhibit some numerical experiments showing the convergence. 
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2. The scheme 

We start by introducing the necessary notation. Derivatives will be approximated 
by finite differences, and the basic quantities are as follows. For any function 
p:t^lRwe set 

D±p(x) = ±^-(p(x±Ax)-p(x)), and D = ± {D+ + D_) 

for some (small) positive number Ax. If we introduce the average 

p{x) = ^ (p(x + Ax) + p(x - Ax)) , 

we find the Leibniz rule as 

D(pq) = pDq + qDp, 

D±(pq) = S ± pD±q + qD±p = S ± qD±p + pD±q. 

Here we have defined the shift operator 

S ± p(x) — pix ± Ax). 

We discretize the real axis using Ax and set Xj = jAx for j £ Z. For a given 
function p we define pj = p(xj). We will consider functions in £ 2 with the usual 
inner product and norm 

{p,q) = Ax^Pjlj, \\P\\ = (P,P) 1/2 , P,qef. 

In the periodic case with period J the sum over Z is replaced by a finite sum 
j = 0, . . . , J — 1. Observe that 

NL = ^b J -|< x ^lbl|. 

The various difference operators enjoy the following properties: 

(P, D ± q) = - (D T p, q) , (p,Dq) = -(Dp,q), p,q£i 2 . 
Further useful properties include 

(u,D 2 + D_u) = ^(u,D_D 2 + u) - ^(u,D + D 2 _u) 

= l{u,(D_D 2 + -D + D 2 _)u) 
(2.1) 2 

= -(u,D_D + (D_-D + )u) 

= ^\\D + D.uf 

since {u 1 D + D 2 _u) = — (u, D_D^u) (because (u, v) = (v,u)) in the first line, and 

D_-D + = -Ax D + D_ = -Ax £>_£>+. 

We also need to discretize in the time direction. Introduce (a small) time step 
At > 0, and use the notation 

for any function p: [0,T] -> R. Write t n = nAt for n £ N = N U {0}. A fully 
discrete grid function is a function u& x : At No — > ffi z , and we write u\ x (xj,t n ) = 
m™. (A CFL condition will enforce a relationship between Ax and At, and hence 
we only use Ax in the notation.) 
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We propose the following implicit scheme to generate approximate solutions to 
the KdV equation ( 1.1 1 

(2.2) u] l+1 =u n 3 - Atu^Du]- At£>^D_u™ +1 , n e No, j 6 Z. 

For the initial data we have 



w. 



° = u (xj), j e Z. 



j 

Remark 2.1. T/izs scheme can be reformulated as an operator splitting scheme as 
follows. Set 

n+l/2 1 / „ „ \ At ( 1 



u 3 2 v ^ +i ' ~i- L > 2Ao 



i.e., u™+ 1 / 2 is solution operator of the Lax-Friedrichs scheme for Burgers' equation, 
applied to u n . Then 

u n+l _ u n+l/2 



-DlD-U 



n+1 



At 

i.e., u n+l is the approximate solution operator of a first-order implicit scheme for 
Airy's equation u t +u xxx = 0. If we write these two approximate solution operators 
as S^ t , and S^ t , respectively, the update formula (|2.2[) reads 



u n+1 = {s£ t o si f ) 

The convergence of this type of operator splitting using exact solution operators have 
been shown in [5J [TT] , with severe restrictions on the initial data. The results in this 
paper can be viewed as a convergence result for operator splitting using approximate 
operators with less restrictions on the initial data, but with specified ratios between 
the temporal and spatial discretizations (CFL-like conditions). 

3. Convergence for smooth initial data 

To show that the implicit scheme can be solved with respect to u™ , we proceed 
as follows: Write the scheme as 

(1 + AtD 2 + D_)u] +1 = u] - Atu^Du™ 

and hence 

((1 + AtD 2 + D-)u n+ \u n+1 ) = \\u n+1 f + At{D%D_u n+1 ,u n+x ) 

= \\u n+l \\ 2 + ^AtAx\\D+D-u n+1 \\ 2 
> \\u n+1 \\\ 

thus 

< |[ (1 + AtD\D-)u n+1 \\ = \\u n - Atu n Du n \\. 
The fundamental stability lemma reads as follows. 



Lemma 3.1. LetUj be a solution of the difference scheme (2.2 1. Then the following 
estimate holds 

(3.1) 

\\u n+1 \\ 2 + AtAx 1/2 ^Ax\\\DlD_u n+1 \\ 2 + Ax 1/2 \\D + D_u n+1 f+ j \\Du n \ 



<IKII 2 , 



provided the CFL condition 

(3-2) A|K|| Q + Ia|| u °||) < 5 6(0,1), 



holds where X = At/ Ax 



3/2 
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Proof. For the moment we drop the indices j and n from our notation, and use 
the notation u for u™ where j and n are fixed. We first study the "Burgers" term 
AtuDu. Let u be a grid function and set 



(3.3) 



w = u — At uDu. 



If the timestep At satisfies (3.2 1 then we have the following "cell entropy" inequality 

.Ax 2 



(3.4) 



1 



-w* < —u 2 



^ Du ^ S ^(Duf, *e(0,l). 



To prove this we multiply (3.3) by u to find 

-w 2 = -u 2 — At-uDu 2 A — (to — u) 2 
2 2 2 2 V ; 



= -u 2 - At l uDu 2 + 1 At 2 u 2 (Duf 
2 2 2 v ; 



i(«). 



Now we have that 



and 



\(a + b) (a 2 -b 2 ) = \(a*-b*)-^{a-bY, 



\a + b) 2 - l -(a 2 + b 2 ) = -\{a-b) 2 . 



4 V ' 2 
For a grid function, this implies 



uDu 1 = -Du A 



2Ax z 



u 2 — u 2 



-Ax 2 (Du) 



(DuY 



Therefore 

1 2 

2 W 



1-^ At 



—Vr 
2 



—vr 
2 



Dm 3 + -AiV (Du) 2 
3 2 V ; 

At „ , „Ax 2 ,„ 2 



AtAx 2 



(Dtz) 3 



Ax 2 



-Dm 3 - 5- 
3 2 



(Du) 



Ax 2 (Duf (^Ax 3 ' 2 Du + ^AxAV - - 



1-5 At , 

- 2 3 



+ Ax 2 (D.) 2 Qax 1 / 2 NL + ^AxA 2 



zZ 2 



1-5 At 3 Ax 2 > 
< -tr ttDu — o- 



(Du)' 



Ax 2 (Duf 



A , , 1 , 9 ,i n2 I — S 

3 IMI + 2 A INI y~ 



< l -u 2 -^Du 3 -S^(Duf, 5 6(0,1), 



where we have employed that A < since A satisfies the CFL condition (3.2 1. 
Estimate (3.4) follows. 

Summing (3.4) over j we get 

(3.5) |H| 2 +<5Ax 2 ||Dh|| 2 < ||u|| 2 . 
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Next we study the full difference scheme by adding the "Airy term" At D\D- 
Thus the full difference scheme (2.2 1 can be written 

v = w- At D\D-v. 

Writing this as w = v + At D+D-v, we square it and sum over j to get 

|H| 2 = |H| 2 + 2At(v, D 2 + D-v) + At 2 ||L> 2 L>_?;|| 2 
= |H| 2 + AtAx \\D_D + v\\ 2 + At 2 \\D 2 + D^v\\ 2 , 



.n+l 



(3.6) 



using the identity (2.1). 

For the function u n this means that 



w\ 



.n+l I 



(3.7) 



+ AtAa; 1/2 (AxA \\D 2 + D_u n+1 f + Ax 1/2 \\D+D-u n+1 \ 
< ||u' l || 2 -<5Aa; 2 ||L»u|| 2 , 



using (3.5 1. This implies 

,2 



(3.8) ||w" +1 || + AtAx 1/2 (Aa;A \\D 2 + D_u n+1 \[ 

+ Ax 1/2 \\D + D- 



-ill 2 ^ 



\Du n f) <iKir 



□ 



Next, we consider what corresponds to the finite difference scheme satisfied by 
the time derivative of the original scheme. 



Lemma 3.2. Letu™ be a solution of the difference scheme (2.2 1. Then the following 
estimate holds 



I D\ u n 1 1 2 + At 2 1 1 D 2 , D_ D\u n 1 1 2 



(3.9) 



AtAx \\D>D-D\u n \\ 2 + SAx 2 WDDiu 11 ' 1 ^ 



provided At is chosen such that 



(3.10) 

Proof. Introduce 



6|| MQ || 2 A 2 + |M|A< 



1-5 



2 ' 



>- 1 \\*(l + 3At\\Du n \\ O0 ), 
5 6(0,1). 



a n = D l ,u n - 1 = - u n ~ v ), neN. 

At y 



Using (2.2 1 we see that this grid function satisfies 

a n+1 = a n - At (a n Du n + u n ~ l D a n ) + At 2 a n Da n - AtD 2 + D_a 
(3.11) = a" - At D{u n a n ) + At 2 a n Da n - AtD 2 + D_a n+1 , n e N. 
Introduce 



(3.12) 



Ar 

P = a-At D(ua) + —Da 2 , 



which means that (3.11 ) can be written as 
(3.13) 



n+l _ 



13 - AtDiD-a 



n+l 



We proceed as before and square (3.12) to find 
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Ax 2 

- At (uaDa + a 2 Du) + At 2 a 2 Da — (Da) . 

We have that 

l - (l)(ua) - ^Da 2 ^j < (D(ua)) 2 + At 2 a 2 (Da) 2 

< 2u 2 (Da) 2 + 2a 2 (Duf + At 2 a 2 (Da) 2 , 

1 Ax 2 

a 2 Da = -Da 3 (Da) 3 , 

3 3 V ; ' 

1 1 Ax 2 

uaDa + a 2 Du = -D (ua 2 ) + -a 2 Du — (Da) 2 Du. 



Using this 



1 „o l^r At „ / o 2At ,\ Ai ,„ AtAx 2 



2 2 < ^a 2 - —D ( ua' - ~—a 3 ) - — a z Du + (Da) 2 Du 



(3.14) + At 2 ( 277- [Da r + -irriDur + A/-V ! Z)n r - ^ ■ 1 

Ax2 fr> ^ 2 
-— (Da) . 

Now we must balance the positive terms with Ax 2 (Da) 2 . To this end we estimate 
Ax 2 (Da) 2 < a 5 , 
At 2 a 2 (Du) 2 < AAa; 1/2 At||ti n || 00 c?|IHt n |, 
At 2 a 2 <2||u n ||' +2||u n - 1 f ) 



OC 



Using these in (3.14| we find 

1„9 l-o Ai / 9 2At o\ .At At,-=-, . A 2 v 
^ 2 < 2 a 2 -y5 N --g-« J +( T + T )« 2 |^|+2At 2 a 2 (^ U ) 



+ A 2 Ax 3 (Da) 2 (V + 2 \\u n \t +2 H^" 1 ^ + (|K 
Ax 2 .„ . 2 



A* / , 2At 



< ia 2 - [ua 2 - ^pa 3 j + At (l + AAx 1 / 2 \\u n \\Jj a* \Du n \ 

+ X 2 Ax 2 (Da) 2 [2Ax (2 \\u n \\l + Wu^Q + ^ (H^U^ + IK" 1 ) 



< -a 2 - —D (ua 2 - —aA + At(l + A ||u ||)a 2 \Du n \ 
2 2 V 3 J 

+ A 2 Ax 2 (Da) 2 [% || U0 || 2 + ^ Huoll - ~S^f (Da) 2 

< I^2_^D ( ua ^ 2 ^L ( A + A t 3 —^a^\Du n \ 
~ 2 2 V 3 J 2 1 1 
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+ Ax z (Da) z I 6||ii || 2 A 2 + ||u || A 



l-8\ ~Ax 



S—(Da) z , 5 6(0,1). 



Here we have enforced the CFL condition (3.10), which in particular implies that 
1 1 t*o 1 1 A < (1 — S)/2, To simplify the numerical expressions, we have employed | < 1. 
Now we multiply with Ax and sum over j to obtain 

(3-15) I \\/3f + ~S^f \\Daf <h\a\\ 2 + ^At WDu^ \\a\\ 2 . 



Writing equation (3.13) as 



1 = (a n+1 + AtD 2 + D_a n+1 y 



we find 



11/311 



\a n+1 \\ 2 + 2At(a n+1 ,D 2 + D^a n+1 ) + At 2 \\D 2 + D_a n+1 f 
|a n+1 || 2 + AtAx ||L> + L>_a n+1 || 2 + At 2 \\D 2 + D^a n+1 \\ 2 . 



Combining this with (3.15) we find 



(3.16) ||a n+i || 2 + AtAx ||L> + £>_a n+1 || 2 + At 1 \\D'j_D. n 



"t'li - ... \/\,.Nd /) „»- Ml- _i_ \+ 2 ll£>2£)_ a »»+l|| 2 

2 ^ /-i i o a i ll n_ nl \ 1 1 _ n 1 1 2 



+ SAx' \\Da n f < (1 + 3At ||£>u n |loo) \\a r > 



□ 



At this point we recall the inequality (cf. Lemma [A. 1[ ): 
(3.17) ll-D«IL < e ||D+i?-u|| + C(e) ||u|| , 

where e is any constant, and C(e) is another constant depending on e. 
The definition of u n , (2.2), can be rewritten 



(3.18) 



y n+l = D t n 



1 

2\x 



D+D_u n - u n Du n - Dj_D-u 



n+l 



where /i = At /Ax 2 — X/Ax 1 ^ 2 . Therefore (using Lemma 
mate) 

1 



3.1 



in the second esti- 



DiD_u n+1 \\ ■ \\a 



v n+l I 



\u n Du ni 



'+ 



2fj, 



D,D_u r - 



1 



< ||a n+1 || + WDu^ \\uo\\ + — {s ll^-^H + C(e) \\u Q \\) 

< \\a n+1 \\ + \\u \\ ( £l \\D 2 + D_u n \\ + C( £l ) \\u \\) 



2^ 



-e D±D-U r ' 



1 

2fj, 



C(e) \\u \\ 



' Wei 



n+l I 



£ i ll u o|| 



Ax 1 ' 2 
2A 



-e \\Di_D-u" 



a 



+ C( £l )|| Uo || 2 + ^C( £ )||« || 

S v ' 

n+1| | + - ||Z?:j__D_u n || + A (choosing E\ and e such that this holds) 



.n+l I 



\DiD- (u n+1 -Ata 



A 



a 



n+l I 



D%D_u 



n+l I 



1 



At\\D±D. 



v«+l I 



A 
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< ||a n+1 || + ^ \\D\D-u n+1 \\ + i ||a"|| (1 + 3At WDu"^) 1 ' 2 + A 

< \\a n+1 \\ + l - \D\D_u n + x \ + i ||a n || (1 + 3A|| Mo ||) 1/2 + A, 



where we have used (3.16) to estimate At ||D^.D_a n+1 1| . Hence 

(3.19) \\D 2 + D-u n+1 \\ <c + ci||a ,l+1 ||+c 2 ||a n ||, 

for some constants Co, c\ and c 2 that are independent of Ax. Exploiting this and 
the inequality (cf. Lemma [A. 1[ ) in (3.16), we get 

|| a »+i|| 2 < || a n||2 + At ( £ \\D\D_u n \\+C(e) \\u n \\) \\a n f 

< \\a n \\ 2 + CAt (e (c + a \\a n \\ + c 2 {{a^W) + (7(e) Kll) ll«™f ■ 

Since is bounded by ||tio||, 

(3.20) ||a" +1 H 2 < ||a n || 2 + Ai(d 1 || a ' l || 2 + d2(||a"|| 3 + ||a n || 2 ||a"- 1 ||)) , 

for constants d\ and d 2 which only depend on ||uo|| and A. Set a n = ||a™|| 2 , so that 

a n +i <a n + At (dia n + d 2 (af/ 2 + a n a n ^A^ . 
Now let A = A(t) be the solution of the differential equation 

^ = d 1 A + 2d 2 A 3 / 2 , A(t 1 )=a 1 >0. 
at 

This solution has blow-up time 

T°° =t 1 + Alnfl 1 dl 



d 1 \ d 2 ^a 0/ 

Furthermore, for t < T°°, A is a convex function of t (since the second derivative 
clearly is non-negative). We now claim that for t n < T°°, we have 

a„ < A(t n ), n e N. 

This holds for n — 1 by construction. Assuming that the claim holds for natural 
numbers up to n, we get 

fln+l < A (t n ) 

+ At(d 1 A(t n ) 



+ d 2 (A(t n f 2 + A(t n )A(t n -i) 1/2 )) 
< A (t n ) + At (d x A (t n ) + 2d 2 A (t n ) 3/2 ^ 



< A(t n+1 ). 
The last inequality follows from 

A(t n+l )-A(t n ) = I A(s)ds 



d 1 A{s)+2d 2 A(s) 3/2 ^j ds 



> jf A(t n ) + At(diA(t n ) + 2d 2 A(t n f /2 ^j ds 

using the monotonicity. Hence, for t < T = T°°/2, ||a™|| < C for some constant 
independent of Ax. 
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Therefore, we can follow Sjoberg |16j to prove convergence of the scheme for 
t < T. We reason as follows: Let u& x {x,t) be the piecewise bilinear continuous 
interpolation 



(3.21) 



MAx(l, t) = vlj + (x - Xj) D+u] + (t - t n ) D\.u] 
+ {x-x j ){t-t n )D\D + v% 



for (x, t) G [xj,Xj+i) x \t n ,t n+ i). Observe that 

UAx(xj,t n ) = Uj, jeZ, neN . 

Note that ua x is continuous everywhere and differentiable almost everywhere. 
The function uax satisfies the bounds 



(3.22) 
(3.23) 
(3.24) 
(3.25) 



\\uAx( ■ ,£)Hz,2(R) < ll w o|| L 2 (R) , 

\\d x u Ax (- ,t)\\ L2{M) < C, 
\\dtu A x{- ,t)\\ L2{M) < C, 



for t < T and for a constant C which is independent of Ax. The first three bounds 
have already been shown, to show the last bound notice that 



£>!£>_ it" < DW 



\Du n \\ < C. 



The inequality (3.251 follows readily from this. 



Then an applica- 
tion of the Arzela-Ascoli theorem using (3.221 shows that the set {uax}ax>o 1S se " 



The bound on dtUAx also implies that ua x € Lip([0, T]\ 



quentially compact in C([0, T];L 2 (R)), such that there exist a sequence {uA Xj } jeN 
which converges uniformly in C([0, T]; L 2 (M)) to some function u. Then we can 
apply the Lax-Wendroff like result from [7] to conclude that u is a weak solution. 



The bounds (3.23), (3.24), and (3.25) means that u is actually a strong solution 
such that (1.1) holds as an L 2 identity. Thus the limit u is the unique solution to 



the KdV equation taking the initial data Uq. 

Summing up, we have proved the following theorem: 

Theorem 3.3. Assume that uq G H 3 (M). Then there exists a finite time T, de- 
pending only on \[Uo\[ H3 ,^\, such that for t < T, the difference approximations 



defined by (2.2) converge uniformly in C(M x [0, T]) to the unique solution of the 



KdV equation \\A) as Ax ^ with At = O (Ax 2 ) . 



Remark 3.4. We can now proceed as in |16) to conclude the existence of a solution 
for all time: We know that the size of the interval of existence [0, T] only depends 
on the H 3 norm of the initial data Uq. But the exact solution of the KdV equation 
preserves this norm, thus we can define the approximations in an interval [T, 2T] , 
starting from the initial value 



1 

Aa; 



lim UAx{x,T)dx, 

Aa;— ^0 



This can be repeated to conclude that there exists a solution for all t > 0. 

Remark 3.5. To keep the presentation fairly short we have only provided details 
in the full line case. However, we note that the same proofs apply mutatis mutan- 
dis also in the periodic case. In particular, the Sobolev estimates provided in the 
appendix are based on summation by parts where the decay at infinity is replaced by 
the periodicity, yielding the same results. 
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4. Convergence with L 2 initial data 



In this section we show that the same difference approximation defined by (2.2) 
converges to a solution of the KdV equation in the case of initial data u £ L 2 (R). 
Clearly we cannot use previous estimates, since those estimates depend on the 
smoothness of initial data. However in [13], Kato showed that the solution of the 
KdV equation possesses an inherent smoothing effect due to its dispersive char- 
acter. In particular, such an effect cannot be present in solutions of hyperbolic 
equations. More precisely, Kato proved that the solution of satisfies the fol- 
lowing inequality: 

/ \u x \ 2 dxdt) <C(T,R), T,R>0, 

which is the main ingredient in the proof of existence of weak solutions of KdV 
equation with initial data uq G L 2 (M). Indeed we prove that the approximate 
solution uaj, lies in 

W = {w e L a (0,T;H\-Q,Q)) \ w t € L 4/3 (0, T; H~ 2 (-Q, Q))} 

which suffices to get compactness in L 2 (0,T; L 2 (—Q,Q)) using the Aubin-Simon 



compactness lemma, Lemma 4.4 



Let the function p be defined as p — p * u>, where 

p(x) = max{l,min{(l + x + R, 1 + 2R)}} , 

and to is a symmetric positive function with integral one and support in [—1, 1]. We 
are interested in this function for arbitrary and large values of R. All derivatives 
of p are bounded. We shall also use that 

r .R 



< -j^p(x) = I lo{x — y)dy < 1. 

J-R 



dx 



Since p is positive we can define the weighted inner product and corresponding 
norms by 

(u,v) p = (u,pv) = Ax^PjUjVj, \\u\\ p = {u,u) p , 

3 

where pj =p(xj). Note that < (1 + 2R) \\u\\ 2 . 

Using summation by parts (recall that (S ± u)j = Uj±\), we have 

(D-D 2 + u,u) = (D_D\u,up) 

= - (D 2 + u,pD + u + S + uD + p) 

= - {D+ (D+u) D+u,p) - (D 2 +u, S+uD+p) 

= -\{D + (D+u) 2 ,p) + ^ ([D 2 + uf ,p) (D 2 +u, S+uD+p) 
= i [(D+uf ,D_p) + ^ ((D 2 + u) 2 ,p) + (D+u,D_ (S+uD+p)) 
= i ({D+u) 2 , D_p) + ^ ((D 2 + u) 2 ,p) + (D+u, DS+uD+p + uD_D, 
(D+uf .--D-P + D+p^j + ~ ({D 2 + u} 2 ,p) + (uD+u,D-D+p) 
(D+uf , \d_ V + D+p^j + ~ ({D 2 + u) 2 ,p) + l - (D+u 2 ,D_D+p) 
A " '(D+u)\D+D_p 
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(D + uf , D_p + ^D+pj + Ax ((D 2 + u) 2 ,p) - l - (u\D + D 2 _p) 
) p = ({D+uf , D_p) + \ ({D+uf , D+p) 



{D_D 2 u,u) = 



So we have 
(4.1) 

Lemma 4.1. Let it!? be a solution of the difference scheme (12. 21). Let N be such 



+ ^\\Dlu\\l- 1 -(u 2 ,D + D 2 _p) 



that NAt = T, and assume that the CFL condition (3.2 1 holds. Then 

N-l 

iKi£+2AtAxx; ( D + u i +1 ) 2 

(4.2) n=0 |jAa;|<K-l 

< |k ||J + 4 M (|k Ar || 2 + || u °|| 2 Uc : 



where /i = At/ Ax 2 = A/Ax 1 / 2 and t/ie constant C depends only on T and uq. In 
particular, for any finite number R, we have that 

N-l 

(4.3) At Ax E {D+u]) 2 <C R , 

n=0 \jAx\<R-l 

where Cr = C{R, ||uq|| , T). 

Remark 4.2. We shall see that this CFL condition is not sufficient to conclude 
convergence of the scheme. For that we need At = O (Ax 2 ) . 

Proof. As before we set 

w = u — AtuDu. 

Set A = At/ Ax 3 / 2 . If the timestep At satisfies the following CFL condition (3.2) 
then we can multiply (3.4) by p to get the "cell entropy" inequality 

(4.4) ^pw 2 <^pv? - ^pDu 3 - S^-p(Du) 2 , 5 6(0,1). 



Summing ( 4.4 ) over j we get 
\x 2 

~2~ H^Hj' ■ 2 " ; . - - ■ 2 



(4-5) \ \\w\\ 2 p + S^f \\Du\\ 2 p < \ \\u\\ 2 p - ^(p, Du*) + (u 2 , D + D.p) 



By (A. 5 ) we have 

\\uD +P \[ x <e\\D_ (uD+p)\\ + C(e) \\uD +P \\ 

< e (\\D + uD +P \\ + \\uD_D + p\\) + C[e) \\uD +P \\ , 

and similarly 

||^-PIL < £ {\\D + uD_p\\ + \\S+uD_D +P \\) + C(e) \\uD_p\\ . 
We use this to estimate 
\(p,Du 3 )\ = \(Dp,u 3 )\ 

< WuDpW^WuW 2 

^liWuD+pW^ + WuD.pWJWuf 

< 2 ( E (\\D+uD_p\\ + \\uD-D + p\\ + \\D+uD + p\\ + \\S + uD_D + p\ 

+ ^(|| Mj D +P || + ||^_p||))|| U || 2 
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< ^s(\\D+uD_p\\ + \\D+uD+p\\)\\u\\ 2 

+ i E (H IP-D+pIL + \\s + u\\ P-d+pIL 

+ ^ ) (hllll^ + PlL + ll«llll^-Plloo))ll"l| 2 

< e (||L> + uL> + p|| 2 + \\D+uD-p\\ 2 ) + A(s, \\u\\) 
<e(((D+u) 2 ,D+p) +((D+u) 2 ,D-p)) +A(e,||« ||) 

where the locally bounded function A now depends on the first and second deriva- 
tives of p. Recall that ||u|| < \\uq\\, cf. (3.1). Hence, 

\\w\\l + 5^f\\Du\\ 2 p <\\u\\ 2 p + A(e,\\u \\)At 

(4.6) + eAt(((D+u) 2 ,D+p) + ((D+u) 2 , D_p)) 

Ax 2 , , \ 
+ — (u 2 ,D+D_p). 

Next we study the full difference scheme by adding the "Airy term" At Z?^Z?_m" +1 . 
Thus the full difference scheme (2.2 1 can be written 

v = w- At D + D_v. 

Writing this as w = v + At D + D^v, we square it, multiply by p and sum over j to 
get 

\\w\\ 2 p = \\v\\ 2 p + 2At (v, D\D_v) p + At 2 \\D 2 + D_v\\ 2 p 
= \\v\\ 2 p + At 2 \\DlD.v\\ 2 p 

+ 2At ((D+v) 2 , £>_p) + At {[D+vf , D+p 
+ AtAx \\D 2 + v\\ 2 p - At (y 2 ,D + D 2 _p 
Combining this with ( |4.6[ ) we get 

\\vf p + At ({D+vf , D+p) + 2At ({D+vf , D. 

+ At 2 \\D 2 + D_v\\l + S^f \\Du\\ 2 p + AtAx \\D 2 + v\\ 2 p 



< \\u\\; + eAt [(D+uy , D_p) + eAt {(D+u)' , D+p 
+ AtA(e,||«o||) + At (v 2 ,D+D 2 _p) . 
Rearranging and dropping some terms "with the right sign" we obtain 
\\vf p + At(2 - e) {{D+vf , D_p) + At(l - e ) ({D+vf , D + p) 
< \\u\\ 2 p + 2At e ({D+uf - {D+vf , Dp) 

+ At (A(e, \\u \\) + \ ||uo|| 2 + \ \\D + D 2 _p\ 



Next, observe that 



[(D+v) 2 , D+p) >AxY, (D+Vj? > 0- 



\j&x\<R-l 



14 



HOLDEN, KOLEY, AND RISEBRO 



Define the locally bounded function B by B(e,z) — A(e, z) + \z 2 + \\D+D?_p\ 
We choose e = 1/2 and recall that v = u n+1 and u = u n . Then we get 



,n+ 



(4.7) 



1 ||%2AtAx £ (£> + < +1 ) : 

|j'Ax|<fl-l 



< \\u n \\l + At((D + u n ) 2 - (D + u n+1 ) 2 ,D P ) +AtB(e,\\u \ 



This is a telescoping sum, and we choose N such that NAt = T to find 



(4.8) 



i;+AtA-E E (p+vy 1 )' 

n=0 |jAx|<fi-l 



< ||^ + A* -(^ 
From this we can easily conclude the proof of the lemma 



u N ) 2 ,Dp) +TB(e,\\u \\). 



a 



Theorem 4.3. Let {uj} be a sequence defined by the numerical scheme ( |2.2[ ) 7 and 
assume that there is a constant K such that At = KAx 2 . Assume furthermore that 
\\ u o\\l 2 (~r) is finite, then there exist constants C\, C^, and C3 such that 



(4.9) 

(4.10) 

(4.11) 



||WAa;|| i ~( 0!T . L 2(_ 



,Q)) 



<Ci, 



\\ u ^x\\l 2 (Q,T-,H 1 ( — Q,Q)) - C2, 
\dtUAx\\ L 4/3( 0t T;H- 2 {-Q,Q)) - ^3> 



where Q = R—1 and u& x is defined by bilinear interpolation from {it™} , cf. (3.21 1. 
Moreover, there exists a sequence of {Axj}J^ 1 with linx, Axj = 0, and a function 
u G L 2 (0,T;L 2 (-Q,Q)) such that 

(4.12) u Ax - ->u strongly in L 2 {0,T;L 2 (~Q 1 Q)), 



as j goes to infinity. The function u is a weak solution of (1.1 1. 



Proof. We first observe that ||wax|| < H^oll so that (4.9) holds. To that end we first 
recall (3.1) which in particular implies that ||ii n+1 || < \\u n \\. Write now 



UAx 



Ax 



-(Wj+i - Wj), (X,t) G [Xj,Xj+x) x [t n ,t n+ x) 



where Wj = uf + (t — t n )D t + u'j. This implies 



/ \u Ax \ 2 dx = E / 
3 Jx - 

= a,e 

j 

2 11 11 2 
= 3 HI 

<IHI 2 
<IKII 2 - 



x j + 1 



A.) 



-(W j+1 - Wj) 



dx 



2 1 
+ 3 

Az 



Oj+i - Wj) 2 + Wj(w j+1 - Wj)) 



E W J- 



The conclusion follows. 

To show (4.10) we calculate that for (x,t) G [xj,Xj + \) x [t n ,t n+ i 

5 x w A x = D+u'> + (t- t n )D l + D + u] 

= a n (t)D+ul + (1 - a n (t)) D+u n +\ 
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where a n (t) = (t — t n )/At E [0, 1). Using this, we find 



||9 a; MAa;|| i 2( 0iT;i 2(_Q i Q)) — / \\d x UAa:( *)|| L 2 (-Q,Q) & 



<2^Ax £ (D 

\jAx\<Q 



At z 



(t - t n f dt 



D + u" 



_J_ 

] At 2 



tn 
t„+l 



(t n+1 -tfdb 



<3^E A - E ((^0 2 + (^? +1 ) 

n |jAx|<Q 

<c R , 



by Lemma 4.1 This, and the fact that ||uAa;(- ^)\\l 2 {~q Q) — ll w o||, proves (4.10). 
Next, observe that in each cell [xj,Xj + i) x [i„,i„+i) 

(4.13) d t u Ax = + (x - x^D+D^u], 

and from the scheme we have 

At 2 

(4.14) n< ,," — ;-) n „» ,," n,," n ; ■>- „»+i 



We claim that for all sufficiently small Aa; (actually for Ax < 1/3): 

(a) For all n E N , 

\\D.Dlu n \\ H _ 3{ _ QiQ) <C\\D + u n \\ LH _ QtQ) . 

(b) For all n E N , 

||£> +J D_u"|| K _ 2( _ aQ) <C||£> + ^|| i2( _ QiQ) . 

(c) The piecewise constant function u^Du" satisfies 

||u£>u|| L 4/3( ,T ,L 2 (-Q ,Q)) - 

for some constant which only depends on Q, T and ito- 
To prove the first part of the claim, let (j> E Hq(—Q, Q) be any test function 

/ (D_D 2 + u n )(j>(x)dx = E D_D 2 + u] (f>{x) dx 

® \jAx\<Q x i 

E D + u " [ J+ D+D_(j){x)dx 

» I ^* s-\ J X A 



|jAx|<Q 



< e \ d + u j\ r +i \^"^)\ dx 

I • A 1/1-1 iC 7 



|j'Ax|<Q 



E \* 

\jAx\<Q 



x j + 1 



u? / |D + D_0(s)-^"(a;)| dx. 



We start by estimating /I, to that end 



Xj + ± -i rXj + i r-x + Ax pZ 

\D + D^(j){x) - <j)"{x)\ dx < -^—2 



Ax 2 



< 



I nXj + i />x-\-Ax pz 



Ax z 



z-Ax 



\<j}'"{9)\ dOdrdzdx 
y/r - x H^'"!!^^) drdzdx 
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< -Ax 3/2 U"'\\ r2( 



Thus 



II<Ax( Y 3Ax|D+u™| 2 ) ( y W" 



1/2/ x 2 \l/2 



|jAx|<Q liAa:|<g 

< 3Ax \\D + u n \\ LH _ Q Q) W"\\ L 2 { _ QtQ) . 
As to /, we calculate 

i= y \ D +^\ / dx 

\jAx\<Q 

< E \d + u]\Va^u"\\ LHxm) 

\jAx\<Q 

b*|<Q b'Ax|<Q 

= IP+^II^C-Q.Q) II^"IIl 2 (-Q,Q) • 

Therefore (a) follows. Claim (b) is proved similarly. 
To prove (c) we first define the cut-off function i] as 



r){x) 



1 \x\<Q, 

|x|>Q + l, 

x + Q + 1 xe[-(Q + l),-Q], 

^Q + l-x xe[Q,Q + l], 



and set rjj = T)(xj). Then we have that 
N-l „,„ N-l 

A*E(Ek«l 2 ) <AtEii^"ii 4 i 3 ( E (on? 

n=0 j n=0 \jAx\<R 

N-l 1/3 AT- 

<(AtYhu n \t) (AtAxY Y ( Du 



V 1 1/3/ ^ ^ 9\2/3 

j , 

n=0 n=0 \jAx\<R 

N ~ 1 1/3/ ^ ,\2/3 



<[AtY\\vu n C) (AtAxY E 

n=0 n=0 |j'A:c|<-H 

<(AtEirt) cf. 

n=0 



by Lemma |4~Tj To proceed we use the inequality 

\\v\L<2(Ax Y v ") 1/4 ( Ax E ( D +"if 



|jAx|<.R \jAx\<R 
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which holds for any grid function v such that Vj = for |j'Ax| > R. This can be 
shown as follows: 

i-i j-i 
v j = E ( w fe+i ~ w fe) = Aa: E ( Vk + z;fc+1 ^ fc+ Ax ~ 

k ——oo k— — oo 

^ 1/2 ^ 1/2 

< (Ai 51 K + Wfe-i) 2 ) (as 51 P+^-) 2 ) 

k —— oo A:— — oo 

<V2(|M| a + ||i;|| a ) 1/a ||D + t,|| 
<2|H|||£> + t;||, 
which implies that 

\\v\\ oo <2\\v\\ 1/2 \\D + v\\ 1 / 2 . 
We shall use this for v = rju n , to that end observe 

Ax fe«") 2 <KH 2 : 
b*Ax|<fl 

Ax ^ (^ + ^u J ") 2 <2Ax ^ ((^ +% ) 2 + ( ) , j+1 D +U «) 2 

|iAs|<fl |jAx|<fl 

<2|| U "|| 2 + 2Ax E ( D + U l) 2 - 

\jAx\<R 

Hence 

\\m n \\t<c\\u n \\ 2 {\\u n f + A X E 

|jAx|<i? 

Thus 

N—l 2 ,g N-l 2 ,g 

a*E( a * E l«wi a ) < A *E( A *Ek^< 2 " " 

n=0 |jAas|<Q «=0 j 



,.\2 



<c R UtJ2\Kf(\Kf + &x E 



|jAx|<_R 



for a constant C depending only on ||ito||, R, and T. This proves (c). 
Now (a), (b) and Lemma 4.1 mean that 



N 

5 2 n|| 4 /3 -A*VMln r> 2 ™|| 4 / 3 

'+ U IIl 4 / 3 (0,T,H-3(-Q,Q)) — ^ r II + UH-^i-Q.Q) 

n=0 



' 2/3 



<CT 1 /3(AtEl|£ , +« n ll^(-o,(3) 



and that 



71=0 

<C, 



_D_Z? + U \\ L 4/3(nru-3(_ n-\\ — 



L 4 / 3 (0,T,ff- 3 (-Q,Q)) 

Similarly, u n Du n G L 4/3 (0, T, L 2 (-Q, Q)) C L 4/3 (0, T, H~ 3 (-Q, Q)). Therefore, 
by gig ), Z^u? G L 4/3 (0,T, J ff- 3 (-Q,Q)). Next, let 



( x ) = Ax ~ ^')X[x 3 -,x 3+1 )(^)- 
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Then (4.131 reads 



d t u Ax = aD t + u n + (1 - a)£» + 5+? 



Therefore 

\\d t UAx\ 



L4/3(0,T,H- 3 (~Q,Q)) - ll a lli~(R) l|- D + U ' l |li*/3(0,T,fl--3(-Q,Q)) 

+ ||1 — a||ioo( 8 ) ||^+^ +U '1lL 4 / 3 (0,T,ff- 3 (-Q,Q)) 

<r 9 II n* ii™ II <r n 

- Z ll ly +" lli*/3(0,r,H-3(-Q,Q)) - °' 



which is ( |4.11 1. 

Using (4.9 1, (4.10), and ( 4.11[ ) we can apply the Aubin-Simon compactness 
lemma (see Lemma |4.4[ ) to conclude that u Ax has a subsequence which conver- 
gences strongly in L*(0,T; L 2 (-Q,Q)), i.e., ( |4.12| holds. 

Note that this is enough to pass to the limit in the nonlinearity. This means we 
can apply the Lax-Wendroff like result of [7] to conclude that the limit is a weak 
solution. □ 

Lemma 4.4 (Aubin-Simon). Let X, B, Y are three Banach spaces such that X C B 
with compact embedding and B C Y with continuous embedding. Let T > and 
{u n }n£N t> e a, sequence such that {it n }„ e N is bounded in L p (0,T; X) and {9 t it n }„ e N 
is bounded in L q (0, T; Y), for any 1 < p, q < oo. Then there exists u € L p (0, T; B) 
such that, up to a subsequence, 



u„ — > u 



L p (0,T;B). 



5. Numerical examples 

We have tested the scheme for two examples where the solution is known explic- 
itly, and for one example where the solution is not known, but the initial data has 
a singularity, and is in L 2 . 

5.1. A one-soliton solution. The KdV equation ( |1.1[ ) has an exact solution given 
by 

(5.1) wi{x,t) = 9(l-tanh 2 U/%/2{x - 3*))) . 

This represents a single bump moving to the right with speed 3. We have tested 
our scheme with initial data Uq(x) = Wi(x,— 1) in order to check how fast this 
scheme converges. In Figure [T] we show the exact solution at t = 2 as well as the 
numerical solution computed using 1000 grid points in the interval [—10,10], i.e., 
Ax — 20/1000. We have also computed numerically the error for a range of Ax, 
where the relative error is defined by 

„ 1nrt T,jLl\ w -i( x jA)-UAx(Xj,2)\ 

a, = iuu jy . 

E^i^iOj' 1 ) 

Recall that we are using Wi(x, —1) as initial data, so that u>i(x, 1) represents the 
solution at t = 2. In Table Q] we show the relative errors as well as the numerical 
convergence rates for this example. The numerical convergence rate indicates that, 
as expected, the scheme is of first order. Note also that we have to use a rather 
small Ax in order to get a reasonably small error. Computing soliton solutions 
is quite hard, since these solutions are close to zero outside a bounded interval, 
and the speed of the soliton is proportional to its height. Therefore, if a numerical 
method (due to, e.g., numerical diffusion) does not have the correct height, it will 
also have a wrong speed. Thus after some time, it will be in the wrong place and 
the error is close to 100%. 
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U 



3 - 



Initial data 

Numerical solution 

- Exact solution 



10 -8 -6 -4 -2 

X 



2 4 6 



Figure 1. Initial data, and the exact and numerical solutions at 
t = 2 with initial data w\(x, — 1) with N = 1000 grid points. 



N 


E 


500 


51.2 


1000 


31.4 


2000 


17.6 


4000 


9.4 


8000 


4.9 


16000 


2.5 



0.70 
0.83 
0.91 
0.95 
0.96 

Table 1. Relative errors for the one-soliton solution. 



5.2. A two-soliton solution. 

soliton, 
(5.2) 



Another exact solution of ( 1.1 ) is the so-called two- 



W2(x, t) = 6(6 — a) 



6csch 



26*)) +asech 2 L 



fa tanh I 



2at)\ - \Acoth U/bj2{x - 26*)^' 



for any real numbers a and b. We have used a — 0.5 and 6=1. This solution 
represents two waves that "collide" at t = and separate for t > 0. For large \t\, 
W2( ■ , t) is close to a sum of two one-solitons at different locations. 

Computationally, this is a much harder problem than the one-soliton solution. 
As initial data we have used uq(x) = W2(x, —10). In Figure[2]we show the exact and 
numerical solutions at t = 20. Although we have used 4000 grid points, the error 
is a staggering 140%! We see that the qualitative features are "right", in the sense 
that the larger soliton has overtaken the slower one, but neither their heights nor 
their positions are correct. For sufficiently small Aa;, the numerical solution will 
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Figure 2. Initial data, and the exact and numerical solutions at 
t = 20 with initial data W2(x, —10) with A" = 4000 grid points. 

be close to the exact also in this case, but it is impractical to calculate numerical 
convergence rates since the computations would take too much time. 

5.3. Initial data in L 2 . We have also tried our scheme on an example where the 
initial data is in L 2 , but not in any Sobolev space with positive index. Furthermore, 
note that all the conclusions in Section 2] remain valid if we restrict ourselves to the 
periodic case. Therefore we have chosen initial data 

!0 x < 0, 

a;" 1 / 3 Q<x<l, 
x > 1, 

if x is in [—5, 5], and extended it periodically outside this interval. In this case we 
have no exact solution available. Therefore we can only determine the convergence 
by viewing solutions with different Ax. In Figure [3] we have plotted the numerical 
solutions at t = 0.5 using 3750, 7500, 15000 and 30000 grid cells in the interval 
[—5,5]. From this figure we can observe that the numerical solutions seem to 
converge nicely to a (smooth) function. The coarser features are already resolved 
using 3750 grid cells, and only the finer structures become more apparent for smaller 
Ax. 

Appendix A. Sobolev inequalities 

For the convenience of the reader we include proofs of the discrete Sobolev in- 
equalities, that are frequently used, but rarely proved. 

Lemma A.l. Given mi, rri2, m^, ni, ni S No, we define m = mi + m^ + m% and 
n = ni + U2- Assume m < n. Consider u G £ 2 (M). Given a positive e. The 




following estimates hold 



(A.l) 
(A.2) 



\\D^D™*D ms u\\ 2 < e \\D^D^u\\ 2 + C(e) \\u\\ 2 
\D™ l D™*D m3 u\? < e WD^D^uWl + C(s) \\u\\ 



2 

2 ' 
2 

2 - 



for some function C(e). 

The same estimates hold in the periodic case where u £ is such that there 

exists a period J € N such that Uj+j = Uj for all j 6 Z, and the norms are taken 
over the period. 

Proof. We here treat the case of the full line only. Assume first that — 0. The 
proof follows by induction. Let m — 1 and n — 2. Then we have 



\d±u\\; 



(u,D + D-u) 
< e \\D+D_u\\l + C(e) \\uf 2 . 



Since \\D + D_u\\ 



\D+u\\ 2 = \\D_u\\ 9 , w e have shown (A.l) in the case with 



m = 1 and n = 2. Assume now that (A.l I holds for all cases with m < N for 
some fixed but arbitrary N, and all n = m + 1. Given mi, mz, ni, n-i such that 
m = mi + mi = N + 1 and n = ni+7i2 = m + 1. We then find 



< e ||r>+ 1 £>" 2 u|i: + C(e) llD™-™ 2 //"-™! 



u 



12 

< e + C(e)(e 1 \\D^D m -u\\ 2 2 + C(e x ) H*), 



22 



HOLDEN, KOLEY, AND RISEBRO 



using the induction hypothesis since m — n-i + m — n\ = to + (m — n) = N and 
m = N + 1. We can rewrite this as 



(1 - C{e)e 1 ) \\D^D^ U \\ 2 2 < e \\D n ^D n Ju\\l + C{e)C{e x ) \\u\\ 2 . 

Given e we choose E\ such that C{e)e\ < \, which proves the case with to — N + 1 
and n = m + 1. By induction we have shown ( A.l ) in all cases where n = m + 1. 

Next we show how to extend this result to n = to+2 for arbitrary to. The general 
case of n > to follows similarly. Let now m, ri2 be such that n = rti + n 2 = to + 2. 
We now have 



\Dl ll D m2 u\\l < e WD^D^uWl + C(e) \\u\\ 2 2 

< e{e x \\D^D^u\\ 2 2 + C( £l ) \\uf 2 ) + C(e) ||u| la 



2 ° ii- 1 -^-^ ' -* H2 V / II 112 

2 , ^ \ 11 i|2\ , r~ir s 11 1 1 2 



= ee x \\D^D^u\\ 2 2 + (eC(e 1 ) + C(e)) \\uf 2 , 

using first that n\ — 1 +712 = n — 1 = to + 1. This proves ( |A.l | in the general case 
with to 3 = 0. 

For an arbitrary TO3 € N we observe that 

^ 3^ jjni! jjm 2 jjm 3 _ i£- 7713 ^ ' / m 3 j ^-j mi +fc^m 2 +m3 — 

fc=0 ^ ^ ' 

using 13 = j(-D_|_ + DJ), which reduces this case to that with TO3 = 0. 
Consider now the inequality (A.2). Observe that 



"j 



j'-i 3-1 

2 _ A„ 



Ax ^ D + U \ = ^ x X! ( u k + Uk+i)D + u k 

k ——oo k— — oo 

which implies that 



l U HOO 



Thus 



I ((u + j < e ||D+«||a + C(e) ||u||a = £ ||Z>±u||a + C(e) |M| 2 



Iz^r^ull 2 < e llD+p^ 1 ^™ 2 ?/)!! 2 + C{e) \\D n ^D^uf 



< eie^D^D^ u^l + Cie^WuWl) 

+ C(e)(e 1 WDl'D^uWl + Cie^WuWl) 
<(e + C(e)) £l \\mpD?u\\l + (e + CT(e))C(ei) |M| 2 . 

(In the rare case that n\ = mi + I and n 2 — m 2 we do not change the first term.) 
Given e we choose E\ such that {e + C{e))e\ < e. This completes the proof of (A.2 1. 

The proof of (A.2) requires some modifications in the periodic case. Let to be 
such that \u m \ = minj=o,...,j-i \v,j\. For j > to we have 

j-i j-i 
" 2 = u 2 n + Air -D+m 2 = u 2 „ + Ax^ ( u k + u k+l )D+u k . 

k—rn k—m 

Thus 

<min|uj| + j ((u + -D+w) | 



I |2 ^ -i |2 

max M 9 - ■ " 



j 3 

' 11 2 , 11 ,-, 1 1 2 , Pi, s 11 1 1 2 
1 1 7/1 - - 

L 



< T ||u|| 2 + £ ||£> +U ||^ + C( £ )|| U || 2 
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Here we have used that 



<e\\D ± u\\l+C{e)\\u\\l 



j-i 



JAxmin \uj\ 2 < Ax \uk\ 2 — \\ u \\ 2 • 

j z — / 

fe=0 
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which implies 



minimi 2 < - ||u||2, 



where L — J Ax remains finite and nonzero as Air — > 0. □ 



Remark A. 2. We could equally well have written the inequalities (A.ll, (A.2| as 

(A.4) WD™ 1 D m2 D rn Hi\\ 2 < e \D\ l D n ^u\ 2 + C(e) \\u\\ 2 , 

(A.5) ||D™ 1 £>" l2 _D m3 u|| < e||L»^ 1 L'!! 2 u|| + C(e) ||u|| 2 . 
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